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Abstract 

We discuss the quantum transport of electrons tlirougli a resonant tunnel junction coupled to 
a nanomechanical oscillator at zero temperature. By using the Green's function technique we 
calculate the transport properties of electrons through a single dot strongly coupled to a single 
oscillator. We consider a finite chemical potential difference between the right and left leads. In 
addition to the main resonant peak of electrons on the dot, we find satellite peaks due to the 
creation of phonons. These satellite peaks become sharper and more significant with increasing 
coupling strength between the electrons and the oscillator. We also consider the energy transferred 
from the electrons to the oscillator. 
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I. INTRODUCTION 



In recent years, there has been great interest in quantum transport through both single{l| 
electronic levels in quantum dots and single molecules[2|. Very interesting and significant 
signatures of the electron phonon interaction have been observed experimentallyjs] in cases 
where the electron-phonon coupling is beyond the validity of simple perturbation theory. 
Single molecular electronics [4] has also attracted much attention, motivated both by the 
scientific challenges and by their potential applications in nanoelectromechanical systems 
(NEMS) B, fl. The most interesting challenges are to verify the uncertainty principlefl, 



50,3 



to study the quantum transport properties as atomic scale quantum effects [8|, l9|l become 
more important and significant, and the fabrication of devices on a nano scale that are 
expected to be faster, more reliable and more sophisticated than existing technology. An 



interesting aspect of NEMS is the interplay between electrons, 
of leads to the system[9|, [101]. Inelastic tunneling spectroscopy 



phonons and the coupling 
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12| is the most direct 



way of observing the inelastic effects in the current-vo^ 
NEMS can be used as ultra-sensitive detectors for mass 



tag e (I-V) characteristics of NEMS. 



131 ]. extremely weak forcesjlj], ultra 



small displacements 15|] and even in chemical 161] and biological 17| applications. At very 
low bias voltage resonances occur with the frequency of the nanomechanical oscillator. The 



characteristic frequency associated with these systems is of the order of 1 GHz[18|. Such a 
high resonance frequency is sufficient to enable the cooling [19] of a nanomechanical resonator 
to its ground state: a necessary condition for these measurements and something on which 
experimental effort is now under way. Moreover, the quantum transport requires very highly 
sensitive measurements such as are achievable using single electron transistors (SET) or 
superconducting SETs (SSET). 

Most of the theoretical work on transport in NEMS has been done within the scattering 
theory approach (Landauer) but it disregards the contacts and their effects on the scattering 
channel as well as effect of electrons and phonons on each other [2^. Very recently, the 
non-equilibrium Green's function (NEGF) approach [2 1| has been growing in importance 
in the quantum transport of nanomechanical systems 2^ . An advantage of this method 
is that it treats the infinitely extended reservoirs in an exact way, which may lead to a 
better understanding of the essential features of NEMS. The pioneering work on NEGF 
was done by Datta and Lake 231]. NEGF has been applied in the study of shot noise in 
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chain models 



24| and disordered junctions[25| while noise in Coulomb blockaded Josephson 



junctions has been discussed within a phase correlation theory approach [2 6;]. The case 
of an inelastic resonant tunneling structure, in which strong electron-phonon coupling is 
often considered, a very strong source-drain voltage is expected for which coherent electron 
transport in molecular devices has been considered by some workers 27]] within the scattering 
theory approach. Inelastic effects on the transport properties have been studied in connection 
with NEMS and substantial work on this issue has been done, again within the scattering 
theory approach [2^. Phonon assisted tunneling of non-resonant systems has mostly been 
shown by experiments on inelastic tunneling spectroscopy (ITS). With the advancement 
of modern technology, as compared to ITS, scanning tunneling spectroscopy (STS) and 
scanning tunneling microscopy (STM) have proved more valuable tools for the investigation 
and characterization of molecular systems [28l| in the conduction regime. Recently, phonon 
assisted resonant tunneling conductance has been discussed within the NEGF technique at 



zero temperature 



29]. 



In the present work, we employ the Green's function method in order to discuss the 
transport properties of NEMS. This is a fully quantum mechanical formulation whose ba- 
sic approximations are very transparent, as the technique has already been used to study 
transport properties in a wide range of systems. The main differences between existing work 
and ours are: in most of the existing literature a very large chemical potential difference 
is considered while we are able to include a range from very small to very large. In our 
calculation the inclusion of the oscillator is not perturbative as the STS experiments ^, [l^ 
are beyond the range of perturbation theory. Hence, an approach is required beyond the 
quantum master equation 
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30| or linear response. In addition, we aim in future to apply 



the same methodology to describe more clearly non-perturbative systems such as a quantum 
shuttle 22|. Hence, our work provides an exact analytical solution to the current- volt age, 



conductance, coupling of leads with the system, very small chemical potential difference 
and includes both the right and left Fermi level response regimes. However, our theory 
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331] where the coupling 



does rely on the commonly used wide-band approximation [3ll, 
between leads and dot is taken to be independent of energy. This provides a way to perform 
transient transport calculations from first principles while retaining the essential physics of 
the electronic structure of the dot and the leads. Another advantage of this method is that 
it treats the infinitely extended reservoirs in an exact way in the present system, which 
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may give a better understanding of the essential features of NEMS in a more appropriate 
quantum mechanical picture. 



II. FORMULATION 

We consider a single quantum dot connected to two identical metallic leads. A single 
oscillator is coupled to the electrons on the dot and an applied gate voltage is used to tune 
the single level of the dot. In the present system, we neglect the spin degree of freedom and 
electron-electron interaction effects and consider the simplest possible model system. In this 
initial work we also neglect the effects of finite electron temperature of the lead reservoirs 
and damping of the oscillator. Our model consists of the individual entities such a single 



quantum dot and the left and rig! 
Hamiltonian of our simple system 



it leads in their ground states at zero temperature. The 
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32| is 



^^dot-ph = [eo + a{b^ + b)] cJcq + ujob'^b , (1) 

where eo is the single energy level of electrons on the dot with cj, Cq the corresponding 
creation and annihilation operators, a is the coupling strength between electrons on the 
dot and an oscillator of frequency uq and 6^ , b are the raising and lowering operators of the 
phonons. The remaining elements of the Hamiltonian are 



-f^leads — efcc|,Cfc (2) 

k=L,R 

-f^leads-dot = -T= ^ Vk (^C^Cq + cjcfc j , (3) 



k=L,R 
1 

k=L,R 

where is the total number of states in each lead. The total Hamiltonian of the system is 
thus 

H = -f^dot-ph + -f^lcads + -f^lcads-dot • (4) 

We write the eigenf unctions of i^dot-ph in a /c-space representation as 

^n{k,xo 7^ 0) = A„exp[-|fc^]if„(A;) exp[-iA;xo] (5) 
^n{k,Xo = 0) = Ar,exp[-ie]HM, (6) 

for the occupied, Xq 7^ and unoccupied, Xq = 0, dot respectively, where xq is the displace- 
ment of the oscillator due to the coupling to the electron and Hn{k) are the usual Hermite 



polynomials. Here we have used the fact that the harmonic oscillator eigenfunctions have 
the same form in both real and Fourier space. 

In order to transform between the representations for the occupied and unoccupied dot 
we require the matrix with elements 

Amn = j K{k, Xo ^ 0)^nik, Xq = 0) dk (7) 

which may be simplified 3J] as 

Amn = , ^ . = / exp (-/c^) Hmik)Hn{k)exp (ikxo) dk 
A/vr2™+"n!m! J ' 

^:^exp (-Ix^) (iia:o)"--L— (8) 

for m < n, where Xq = A/uq, A = o? juj^ and i^J^(x) are the associated Laguerre polyno- 
mials. Note that the integrand is symmetric in m and n but the integral is only valid for 
m < n. Clearly the result for m > is obtained by exchanging m and n in equation (JHl) to 
obtain 

I mm[m~^ (Hx,,)^'-"'^ C-x^) (9) 

Amn-^ maxKn]! I 4^oj U^^oj ^^mKn] U V • W 

III. GREEN'S FUNCTIONS AND QUANTUM TRANSPORT 

In order to calculate analytical results and to discuss the numerical spectrum of the 
quantum transport properties of nanomechanical systems, our focus in this section is to 
derive an analytical relation for the current, and the differential conductance, as a function 
of applied bias. In obtaining these results we use the wide-band approximation where the 
self-energy of the dot due to each lead is considered to be energy independent and is given 
by 

,2 deK 



= riD \Vk\ 



00 E - ex 
2 



-27rinD \Vk\ (10) 



-iir 



K 



where nD is the constant number density of the leads, K = L,R represent the left and right 
leads and Tk is the damping factor (F^, = F/j = F). Similarly 

Sx = ra* = +iir^. (11) 



We solve Dyson's equation using i?dot-ieads as a perturbation. For the more general systems 
we aim to treat in future this is a reasonable small parameter. In the present case, however, 
we can find an exact solution. The retarded and advanced Green's functions on the dot, 
with the phonon states in the representation of the unoccupied dot, may be written as 

G^nn' =E^-^-"^^n'-' (12) 
m 

where gn^'' is the retarded (advanced) Green's function on the occupied dot, 

) = [£; - eo - (n + \)uJo + A ± iP] (13) 

with A = /ujq. 

The lesser self-energy may be written as 

E< (£;) = \TfK{E) . (14) 

where fL{R) are the Fermi distribution functions of the left and right leads, which have 
different chemical potentials under a voltage bias, including a contribution from the state of 
the oscillator. 

For the present case of zero temperature the lesser self-energy may be recast in terms of 
the Heaviside step function 9{x) as 

ll<{E)^\Te{eYK + \uJo-E)5n,o, (15) 

where e-px is the Fermi energy on lead and the Kronecker delta, 5„^0) signifies that the 
oscillator is initially in its ground state, n — Similarly one can calculate the greater 
self-energy as 

s> (E) = -ir[i - . (16) 

The lesser Green's function is related to the density matrix through 

= -2iG<„, , (17) 

Here G^^, is the full lesser Green's function including the dot and the leads. With the help 
of the density matrix formulation, the current formula is 

Tr [pi^ = ii Tr (^G< j) , (18) 
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where / is the current operator. Using this formula, we calculate the current from the 
contact through both barriers and the central region with the oscillator on the dot. The 
general expression for the current is 



/ = ^ I {Tr [(E< - S<) (G^- - C^)] + Tr [((S^ - SI) - (S?, - ) G<]} dE . (19) 

We can obtain the same result by calculating the current from the dot into one of the leads, 
which gives 

/ = |- I {Tr[-S< (G"- - an] - Tr[(S?, - S^,) G<]} dE , (20) 

where the first term in the above expression vanishes for n > 0. The lesser Green's function 
in the presence of the oscillator is given by 

Q< ^ Qrj^<Qa ^^^^1 S< = S< + S< . (21) 

Using the same damping factor in each lead (F^ = F/j = F), the final expression for the 
total current (/) reduces to 

J = ^ 1 Tr (S< - S<) (G'- - G") dE . (22) 

IV. AVERAGE ENERGY 

To calculate the energy transferred from the electrons to the oscillator we return to 
equation fl2Ul) and note that the contributions to the trace with n > may be identified 
with the current due to those electrons which have been scattered inelastically with the 
creation of n phonons. As the lesser self-energy factors in the 1st term are zero for n > 0, 
the inelastic contributions are solely contained in the 2nd term. The first term in equation 
f l2Up does contribute to the total current calculated in equation 0221) but does not contribute 
to the energy of the oscillator. We may therefore use the lesser Green's function (or density 
matrix) to calculate the energy transferred to the oscillator to obtain 

Eph = J2 ^^oYn/ 1 where F„ = |i^ f G<„ dE , (23) 
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From equation (123!) we may write the lesser Green's function in terms of the lesser self-energy 
and the retarded and advanced Green's functions as 
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(24) 



where we note that, as we are working at T = the self-energy terms are only non-zero for 
the zero phonon state. Hence we have 



2 47r 



m,k 
X 



A A* 



E - eo - (m + I) u;o + A + ir 
ir^ (epL + 1^0 ~E)+ iVd (epij + - E) 



dE. 



(25) 



E - eo - (A; + i) LJo + A - ir 

We note, however, that this expression is non-zero even when epL = ^fr and eo < ep due 
to the dot being permanently occupied in these circumstances. To remove this term we 
subtract the contribution when the 2 Fermi levels are equal. This reduces the expression for 

Yn to 

^r^e ^^^^ 



2 47r 



E 



A A* 



E_- eo - (m + I) cuo + A + iF 



(26) 



After integrating the above expression [35|], we arrive at the final result 



^ F^e X — -\ An rn.A-r\ ^-An kA-t 



(h - 



2 47r ^ (k- m)uJo + 2iF 

m,k 



X < In 



epL - eo - muo 



A-iF 



epL - eo - kujQ + A + iF 



In 



epR - ep - mujQ + A - iF 
epi?. - eo - kuo + A + iF 



(27) 



Hence, the average energy transferred to the oscillator may be calculated using equation 

(ESD. 

V. DISCUSSION OF RESULTS 



The differential conductance is shown graphically in figure [T] as a function of 



voltage for different values of coupling strength, using the same parameters as 



29 



app 



31 



ied 



37|, 



: the single energy level of the dot eo = 0.5, the characteristic frequency of the 
oscillator uq = 0.3, the damping factor F = 0.3co'o and the chemical potentials < epL < 1 
and epij = 0. These are chosen to illustrate the physics of such systems rather than to 
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Coiipling Slieiigth(ctJ 



FIG. 1: Differential conductance (dimensionless) as a function of applied voltage epi (in arbitrary 
units) and coupling strength a. Gate voltage eo = 0.5, oscillator frequency luq = 0.3, self-energy 
r = 0.3^0. 

represent a specific implementation. The oscillator induced resonance effects are clearly 
visible in the numerical results. It must be noted that we have obtained these results in 
the regime of strong and zero or weak coupling of the oscillator with the electrons on the 
dot. The coupling between the leads and the dot is considered to be symmetric and we 
assume that the leads have constant density of states. With increasing coupling strength, 
the number of satellite peaks also increases while for zero or weak coupling we find only 
the basic resonance. This confirms the effect of the coupling between the electrons on the 
dot and the single oscillator mode where higher energy electrons are able to drop to the 



dot energy by creation of phonons. We note the similarity of figure [T] to figure 3a of 39[ , 
which refers, however, to the transmission amplitude of an interference device, albeit using 
a similar Hamiltonian. Transport processes involving creation or annihilation of phonons 
are a common feature of NEMS. 

Closer analytical examination of the expression for the differential conductance (!26|) shows 
that the main resonance peaks occur when the applied voltage, ei^^^ is equal to the energy 
eigenvalues of the coupled dot electron and oscillator. The main peak {n = 0) is given by 



the Lorentzian form with its center at the epi = ^o — known as a Breit-Wigner 36|, |40L |41| 



resonance. The satellite peaks due to the emission of phonons can be seen on the positive 
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energy side with epi = e^ — ^ + nujQ where ujq is the characteristic frequency of the oscillator. 

The main or basic resonance peak is the elastic or zero phonon transition. The amplitude 
of the satellite peaks or steps is much smaller than the basic resonance peak. The electrons 
that tunnel onto the dot can only excite the oscillator mode as at zero temperature there are 
no phonons available to be absorbed. Moreover, we have seen that with increasing coupling 
strength, the number and intensity of the satellite peaks increases but their intensity always 
remains much smaller than the main peak. The peaks or steps in the current characteristics 
vanish if the upper electrochemical potential is smaller than the dot energy plus the oscillator 
frequency. 



T — « — « — « — « — I — « — « — ' — « — I — « — « — « — « — I — J — « — « — « — I — « — « — « — « — r 




Gate Voltage (tg) 

FIG. 2: Differential conductance (dimensionless) as a function of gate voltage eq, with applied 
voltage epL = !> oscillator frequency ujq = 0.6, self energy T = O.Swq and coupling strength 
a = 0.4t(;o (dotted line), O.Gwq (light solid line), O.Sljo (bold solid line). 

The differential conductance as a function of gate voltage, eo, is shown in fig. [2] for various 
coupling strengths at T = 0. The main peak at eo = epL + A) corresponds to elastic or zero- 
phonon transition and the satellites peaks are due to emission of phonons corresponding to 
n = 1, 2, 3, 4, . . .. This shows more and more satellites corresponding to every multiple of 

With increasing coupling strength while keeping the temperature zero, we see that the 
energy transferred to the oscillator increases with increasing coupling strength, while the 
amplitude of the satellite peaks is much smaller than the main peak which is shifted toward 
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the right by a factor A. The amphtude of the main peak is also affected: its magnitude 
decreases with increasing coupling strength. 
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FIG. 3: Average Energy transferred to the oscillator as a function of applied voltage epL with gate 
voltage eo = 0.5 (left figure) and as a function of gate voltage eo with epL = 1 (right figure). The 
oscillator frequency ljq = 0.1, selfenergy F = OAloq, and coupling strength a = OAloq (dotted line) 
and a = O.Gi^o (solid line). 



In fig. [3|, we plot the average energy that is being transferred to the oscillator per trans- 
mitted electron as a function of applied bias and gate voltage. At zero temperature the 
oscillator can only gain energy from the electrons. We note that there is some structure 
as individual phonons are excited but there is also a saturation level oc a^. The peak just 
below Co ~ epji in the right hand figure is due to the fact that there is no elastic transmission 
in this regime and all transmitted electrons result in the creation of phonons. Note that 
this is the average energy transferred when the system starts in its ground state and should 
not be confused with the energy transferred after many electrons have interacted with the 
oscillator. We consider T = in this work for simplicity and will consider finite temperature 
effects later. Moreover, the phonon energy of the oscillator on the dot and the level width are 
both typically larger than the experimental temperature. Our theory is in good agreement 



with the growing body of theoretical 29|, 
area. 



31 



32 



36| and experimental 37|, |38|| work in this 
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VI. SUMMARY 



In this work, we analyzed the dynamics of a nanomechanical oscillator coupled to a res- 
onant tunnel junction by using the Green's function approach without treating the electron 
phonon coupling as a perturbation. We have derived an expression for the current and dif- 
ferential conductance and discuss it in detail for different values of the coupling strength. 
Wc have found steps/peaks in the current spectrum as a function of the chemical potential 
difference in addition to the main resonant step, due to the transfer of energy from electrons 
on the dot to the oscillator. We have also studied the effect of gate voltage. We also derive 
an expression for the average energy transferred from the electrons to the oscillator. We 
have shown that the steps grow with increasing coupling strength of electrons on the dot 
and the oscillator. This confirms that the additional satellite peaks or steps in the spectrum 
of numerical results are due to the transfer of energy from the electrons to the oscillator. 
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